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Abstract 

We present a calculation of defect-defect correlation functions in the defect 
turbulence regime of the complex Ginzburg-Landau equation. Our results do 
not agree with the predictions of generic scale invariance. Using the topo- 
logical nature of the defects, we prove that defect-defect correlations cannot 
decay as slowly as predicted by generic scale invariance 
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Spatiotemporal chaos occurs in extended systems with many interacting degrees of free- 
dom []]]. Typically, it appears in nonequilibrium pattern-forming systems slightly above 
their threshold of instability 0. The study of spatiotemporal chaos has been advanced by 
the development of experimental systems which are precisely controlled and have a large 
aspect ratio Such systems have large statistically homogeneous regions relatively free 

from boundary effects. A key question to address is whether such regions can be described in 
terms of hydrodynamic-like theories, focusing on collective behaviors and long-wavelength 
descriptions. We wish to address an aspect of this question by considering the coherent 
structures known as topological defects (or spirals or vortices) in a system that exhibits 
spatiotemporal chaos: the complex Ginzburg-Landau equation |J. This equation describes 
the slowly varying amplitude and phase in an extended system which undergoes a Hopf 
bifurcation to an oscillating and spatially uniform or oscillatory and spatially periodic state. 
The equation exhibits many interesting patterns, but we will restrict our investigation to the 
Benjamin-Feir unstable (or defect turbulent) regime 0, in which topological defects occur 
in the context of spatiotemporal chaos. Other systems, such as Rayleigh-Benard convec- 
tion electrohydrodynamic convection in liquid crystals [|J, capillary ripples ||, cardiac 
tissue JS|, and chemical reactions |S[], can exhibit similar defect-turbulent behavior. In this 
letter, we will examine the defect-defect correlation functions and relate them to the ideas 
of generic scale invariance, which is a theory for describing nonequilibrium systems with 
conservation laws. We will see that our results do not match the predictions of generic scale 



invariance 



|10J. Finally, we will prove that the generic predictions cannot be correct for 
topological defect correlations. 

The complex Ginzburg-Landau equation is given by 

d t A = A - (1 + iC )\A\*A + (1 + z6 x ) |^ + (1 + ib y )^ (1) 

where A is a complex field in two dimensions. This equation can have topological defect 
solutions where A = (both lie A and 1mA are zero) ||[nj. These defects can occur in 



either static arrangements or in dynamic ones called defect turbulence, where defects are 



continuously nucleated and annihilated in pairs and are moving about. We wish to focus 
on the latter case, the Benjamin-Feir turbulent instability regime J7|, which occurs when 
l+b a c < 0. In this region of parameter space, all periodic solutions of the complex Ginzburg- 
Landau equation are unstable. For comparison with the ideas of generic scale invariance 
IIJ, we will focus on the anisotropic case b x ^b y . 

We note that the topological defects come in two varieties. The type of defect depends 
on how the phase of A changes as we go counterclockwise once around the defect. A defect 
with a phase jump of 2tt has a topological charge of +1, while one with a jump of — 2tt 
has a topological charge of —1. This is analogous to the right-handed and left-handed 
single-armed spirals in Rayleigh-Benard convection. We let p+(r) equal the density of +1 
defects and p_(r) equal the density of —1 defects. We can then define a "topological" order 
parameter, p(r) = p + (r) — P-(r), which is just the density of the defects weighted by their 
topological charge. This order parameter is conserved in a system with periodic boundary 
conditions: J v p(r)dr = 0, as defects can only be created or destroyed in ± pairs. We focus 
on the order parameter p(r) as an effective coarse-grained field, which we conjecture can be 
described by some hydrodynamic equation of motion. 

In equilibrium systems, spatial correlations typically decay exponentially. For nonequi- 
librium systems (such as those with an external driving force) the situation can be quite 
different. For nonequilibrium systems with a conservation law and external noise, spatial 
correlation functions can decay algebraically. It has been suggested that this algebraic decay 
is expected to occur for a broad range of conditions, and it has been called "generic scale 
invariance" [|Hj. Some extended deterministic chaotic systems also exhibit algebraic decay 



12| |14[ . In at least one of these examples the chaotic fluctuations appear to play the same 



role as stochastic noise |12| . The complex Ginzburg-Landau equation would seem to satisfy 
the criterion for generic scale invariance. It shows nonequilibrium behavior since it cannot 
be derived from an underlying potential (i. e. it is non-relaxational). In a system with 
periodic boundary conditions, the topological order parameter p(r) is conserved. Finally, 
we conjecture that the chaotic noise in our system plays the role of stochastic noise. 



With this set of conditions, we could have a hydrodynamic equation for the conserved 
order parameter of the form 

d t p(r,t) = r{p(r,t)}+rj(r,t) (2) 

where T is a general conserving; operator on p, such as r V 2 + T^V 2 ) 2 + T 2x d* + T 2y d*. It 
can also contain nonlinear terms (e. g. V ■ [(V 2 p) (V p)\ ) . The stochastic noise term rj is 
determined by: 

( V (r,t)) = (3a) 

(r ] (r,t)? 1 (r , ,t>)} = D5(r-r')5(t-t>), (3b) 

where D must be composed of differential operators for our strictly conserved order pa- 
rameter. This conserving noise term represents the effect of the chaotic fluctuations in the 
complex Ginzburg-Landau equation. There is evidence from the mapping of the Kuramoto- 



Shivashinskii equation to the Kardar-Parisi-Zhang equation |15| and from coupled map lat 



tices [0 that this identification of spatiotemporal chaotic fluctuations with stochastic noise 
is not unreasonable. 

For systems with nonconserving noise {i.e. D is a constant), equation (Q) is expected to 
always give rise to power law decays in the two point correlation function G p {r) = (p(r)p(O)), 
as well as in higher order correlation functions. For systems with conserving noise (e. g. 
D = DiV 2 ) the situation is somewhat more complicated [fLOfl . If the system is isotropic, 
then one obtains exponential decays in G p (r), but power law decays occur in higher order 
correlation functions. Systems which break the isotropy give rise to algebraic decay in G p (r). 
For systems with cubic symmetry, one expects G p (r) ~ l/r d+2 for large r. For systems which 
break the cubic symmetry one expects G p (r) ~ l/r d . We will be working in the last regime, 
where in a two-dimensional system with broken square symmetry, generic scale invariance 
predicts that 

G9 enerio {r) = (p(r)p(0) ) „ 1/^ (4) 

4 



for large r ||10|| . We shall compare this prediction to the results from our numerics. 

Before we consider the numerics, we should provide some caveats. The ideas behind 
generic scale invariance depend upon showing that the nonlinearities are irrelevant in a 
renormalization group sense. It has proven notoriously difficult to treat topological defects 
in a perturbative manner. An example of this is the Kosterlitz-Thouless transition MB 



We also note that generic scale invariance requires short-ranged interactions. There is some 



evidence that this is true for the defects in the complex Ginzburg-Landau equation |T7 
but collective effects might be important. Finally, the mapping of chaotic fluctuations to 
stochastic noise could break down. 

We have numerically solved equation ([[]) with periodic boundary conditions in the tur- 
bulent regime using a pseudo-spectral code. The system is 240 x 240 in real space, with 360 
Fourier harmonics in both directions. We use the parameter values c = 1.5, b x = —0.75, and 
by = —3.0. The time step used was At = 0.02. We initially "equilibrate" from a state with 
two oppositely charged defects to a state with fluctuations about some average number of 
defects. This takes typically 5000 time steps. We note that the defects do not form bound 
pairs. When a pair is created, the defects tend to move apart, and when they eventually 
annihilate, they usually do so with a defect other than their initial partner. This illustrates 
that we aren't in a Kosterlitz-Thouless bound pair phase [|l(| . In figure [l] we show a snap- 
shot of part of our system. To find the defects in our system, we examine the change in the 
phase of A as we go counterclockwise around each plaquette on our lattice in real space. 
A change of signifies that the plaquette does not contain a defect, while changes of ±2ir 
reveal that a defect exists in the plaquette. Once we have found the defects, we can then 
proceed to calculate n(t), the total number of defects in the system, and G p (r). To do the 
averaging, we have run for 750, 000 time steps (one month of CPU time on an IBM RS/6000 
model 550: locating the defects is the time consuming part). We only sample G p (r) and 
n(t) every 10 time steps, because adjacent time steps are not statistically independent. We 
have calculated that (n(t)ri(O)) — (n) 2 ~ e - */ T , with r ~ 115 time steps. It has been pre- 
dicted 0] that the probability of finding a particular value of n in the system is given by 



P(n) ~ e ( n ^ 2 / 2 ^ . We have calculated the various moments of our distribution P(n), 
and we find (n) = 422.8±0.3, a 2 = 397±30, and a skewness of 0.014 and kurtosis of -0.026, 
which is in good agreement with the predictions. In figure Q we present the results for G p (r) 
with r in both the x and y directions. For both directions the typical nearest neighbor is 
of the opposite sign: the charges are thus screened. Similar behavior for vortices in random 
wave fields has been observed PDfl . In figure [3] we show log-log plots of |G p (r)|. We also 



show lines that represent the slope |G p (r)| should have if it decayed like 1/r 2 . We note that 
at the right edge of the figure, we have reached the point where our data is dominated by 
statistical noise. It is clear that neither direction shows the expected 1/r 2 decay. Our results 
are at variance with the predictions of generic scale invariance. We expect that the theory 
is not applicable to these sorts of systems with strong constraints placed on them due to the 
topological nature of the order parameter. We can't explicitly determine an equation like 
equation |] for the defects, but we can discuss the results for |G P (V)| from the viewpoint of 
topological constraints. We will show how these constraints place bounds upon the decay 
rate of the correlation function. 

We define the excess order parameter in a region to be 



5 pi 



/ (p(r + r )-, o)d r 

JreB(L) 



(5) 



/rgB(L) 

where B(L) represents a circle of radius L about a point r , which we take as r = due to 
translational invariance. For nontopological objects, the constraint is given by 8pi < a\L 2 , 
where a\ is some numerical constant. The excess of a nontopological object in a particular 
region must scale as the volume of that region, since each individual object occupies a fixed 
volume. For topological objects, this constraint is different, i. e. 5pi < a 2 L, where a 2 
is again some numerical constant. The constraint arises from the fact that any excess of 
topological defects in a region must be detectable simply by traversing the perimeter of 
that region. Each topological defect has an "arm" with characteristic width that must pass 
through the perimeter of the region. Examples of this are the spiral arms of the defects 
in Rayleigh-Benard convection, extra rows of atoms for dislocations in crystals, and in our 
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case lines of Re A = and 1mA = 0. When a region contains the maximum excess number 
of defects allowed, each of these lines takes up a fixed amount of the perimeter of the region. 
Since the maximum excess number of topological objects scales linearly with the number of 
lines, and the number of lines scales as the perimeter of the region, we must have that the 
maximum excess number of topological objects scales as the linear size L of the region. 

If we assume that the correlation function G p (r) decays at the same asymptotic rate 
independent of the direction of r, i.e. G p (r) ~ f(8)g(r), where g(r) ~ l/r a for large r, then 
with this constraint we can show for two dimensions that a must be greater than 2. This 
result also requires that J 2w f(6)d6 ^ 0, which we expect to be true except for special cases. 
An example of such a correlation function satisfying both assumptions is given in reference 
20fl . To show that a > 2, we begin by squaring equation (|j) and averaging the result (over 



the noise or over time and space). This gives us the constraint equation 

(5p 2 L ) = f f drdr 'GJr - r ') < a 2 L 2 . (6) 

Jr£B{L) Jr 'GB(L) 

Separating G p (r — r ') into radial and angular components yields 

Wl) = 7T / f^ d9 / 9(R)w(R)dR (7) 

27T JO Jo 

where we can derive that 

w(R) = AnRL 



(8) 



This result is obtained by considering how often \r — r'\ = R when r G B(L) and r ' G B(L), 
i. e. 

w(R)= f I drdr 'S(R)8(t — r ' — R), (9) 

JreB(L) Jr 'eB(L) 

where S(R) = 2wR (the circumference of the circle of radius R). From this definition 
w(R) can be calculated by writing the 5-function in integral form, and then performing 
the resulting integrals. By assuming that G P (R) exhibits its asymptotic behavior outside of 
some R > r min , we can split the radial integral in equation (0) into two parts: 



/ g(R)w(R)dR = / + / g(R)w(R)dR. (10) 

JO Jr min JO 

We can then examine the large L limit of (Sp 2 L ) for various values of the power law exponent 
a. By rescaling all lengths by 2L and then expanding w(R) about small R, we can show 
that the leading order behaviour of (8p 2 L ) is 

a < 2 : (5p 2 L ) ~ L 2 L 2 ~ a 
a = 2:(5p 2 L )~L 2 \og(L) 

a > 2 : (5p 2 L ) ~ L 2 . (11) 

Within our assumptions, this result means that for topological objects the results of 
generic scale invariance cannot hold, since the constraint given by equation (0) would be 
violated. In fact, what we have provided is a bound on a. For topological objects a must 
be strictly greater than 2. Generic scale invariance predicts a = 2. The simple geometric 
nature of topological objects prevents them from having correlation functions that decay 
as certain power laws. An added conclusion is that if the topological objects form ordered 
states, they must be of the antiferromagentic variety (e.g. alternating + and — vortices) in 
at least one direction, in order to satisfy the topological constraint. An example of such a 
state is given in |pl| . 

In our analysis we have only considered the largest possible fluctuations. We expect 
these fluctuations to be rare, and hence expect a faster decay than the bound we provide. 
As an analogy, for nontopological objects the analysis presented here would predict that the 
correlation function can be at most a constant for large r; in practice systems like spins or 
atoms have connected correlation functions that decay to zero, either as power laws or as 
exponentials. We expect similarly for real topological objects the decay of the correlation 
functions will be faster than the bound given here. Numerically, we have found that we are 
indeed far from saturating this bound. 
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FIGURES 

FIG. 1. Snapshot of a 70 x 70 region. The solid lines are where Re A = 0, and the dashed lines 
are where ImA = 0. Filled circles (•) are vortices with topological charge +1, and the open circles 
(o) have charge —1. 

FIG. 2. G p (r) versus r. The solid line is for the x direction, while the dashed line is for the 
y direction. Also shown is a line for G = 0. Note that G attains its asymptotic limit of from 
different sides of this line. 

FIG. 3. Log-log plot of |G p (r)| versus r. The solid line corresponds to the x direction and 
the dashed line to the y direction. Also shown is a line with slope that would correspond to 
|G p (r)| ~ 1/r 2 . Note also the break in the line for the y direction which corresponds to the zero 
crossing of G. 
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